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We suggest to excite dipole-forbidden transitions in quantum-mechanical systems by using ap- 
propriately designed optical nanoantennas. The antennas are tailored such that their near-field 
contains sufficiently strong contributions of higher-order multipole moments. The strengths of these 
moments exceed their free space analogs by several orders of magnitude. The impact of such exci- 
tation enhancement is exemplarily investigated by studying the dynamics of a three-level system. 
It decays upon excitation by an electric quadrupole transition via two electric dipole transitions. 
Since one dipole transition is assumed to be radiative, the enhancement of this emission serves as 
a figure of merit. Such self-consistent treatment of excitation, emission, and internal dynamics as 
developed in this contribution is the key to predict any observable quantity. The suggested scheme 
may represent a blueprint for future experiments and will find many obvious spectroscopic and 
sensing applications. 
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I. INTRODUCTION 

Optical nanoantennas have changed our perception of 
how light can interact with matter. If such a nanoan- 
tenna is made from a noble metal, it supports local- 
ized surface plasmon polaritons (LSPP) at distinct fre- 
quencies in the visible and near-infrared spectral range. 
LSPPs are excited when the electromagnetic radiation 
is resonantly coupled to the charge density oscillation in 
the metal. LSPPs allow the focusing of light into vol- 
umes (hot spots) inaccessible by classical optical devices. 
Moreover, the hot spot intensities may exceed the inten- 
sity of the external illumination by orders of magnitude. 
These remarkable features render optical nanoantennas 
as prime candidates for controlling and improving the 
interaction of far-field light with other nanoscopic build- 
ing blocks such as quantum dots, atoms, or molecules^"— . 
Here, we shall treat such nanoscopic building blocks at 
a quantum-mechanical level to fully grasp their internal 
dynamics and to appropriately describe their properties. 

Many applications have been developed that exploit 
light interaction with such hybrid systems consisting of 
a nanoantenna and a nanoscopic building block. Piv- 
otal examples can be found in the field of biosensing^ 
and photovoltaicsiSiii. In studies that considered indi- 
vidual hybrid systems, the most important finding con- 
sisted in showing that optical nanoantennas can modify 
the radiative decay rate of emitters^^. By taking advan- 
tage of this effect, one can either enhance the emission 
rate or facilitate non-radiative decays. The interaction 
of light with such hybrid systems was first studied for 
highly symmetric optical nanoantennas as, e.g., metallic 
nanospheres. However, more complex antennas can be 
equally considerediS,. In such cases, the emission charac- 
teristics of these hybrid systems might be entirely gov- 
erned by the optical nanoanantenna and may strongly 
deviate from the radiation pattern of emitters in free 
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FIG. 1. Left: a plane wave (purple) scatters two close sil- 
ver nanospheres. Middle: molecules in-between undergo a 
quadrupolar excitation resulting in a luminescence signal at a 
lower wavelength (green). Right: excitation and emission of 
a single molecule. 



space ^ . These new engineering possibilities pave the 
way for the development of highly directed single pho- 
ton sources and other applications i^ii^ 

However, in all the aforementioned studies the in- 
teraction of optical nanoantennas with the quantum- 
mechanical system has been discussed by resorting to 
more or less appropriate approximations. To be specific, 
apart from a few exceptions^^"-*^*, the effect of higher or- 
der multipole fields in the vicinity of the antenna has 
been restricted to a pure electric dipole field; although 
experimental indications for the influence of higher order 
multipole fields have been reported in the literature i^^i^° 
The restriction to electric dipole fields seemed to be an 
obvious and quite reasonable approximation, since in free 
space the interaction of a quantum-mechanical system 
with higher-order multipole fields is orders of magnitude 
weaker than compared to that of an electric dipole field; 
as discussed in detail below. Moreover, none of the afore- 
mentioned studies considers the modification of the dy- 
namics of the quantum-mechanical system subject to a 
higher-order excitation. However, subsequent calcula- 
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tions will show that this is of crucial importance. 

With the self-imposed limitation to consider the field 
close to an optical nanoantenna as to be electric dipo- 
lar, the potential of optical nanoantennas is unnecessar- 
ily restricted because there is a much larger flexibility in 
tailoring their near-field. In this work the potential of op- 
tical nanoantennas to significantly enhance higher-order 
multipole fields in the vicinity of the nanoantenna is con- 
sidered. The enhancement is such strong that it can be 
essential in the interaction with a quantum-mechanical 
system. Most notably, it enables the significant excita- 
tion of transitions which are typically considered forbid- 
den in free-space. 

To analyse such processes, a theoretical framework nec- 
essary to calculate non-dipolar transitions of quantum- 
mechanical multilevel systems in the vicinity of optical 
nanoantennas is established in this contribution. These 
non-dipolar transitions, often termed as forbidden ones, 
can be tremendously enhanced due to large higher-order 
field componentsi^iiSii^i^. In passing we note that 
these higher-order components may also be used to ex- 
cite dark modes in plasmonic systems using strong field 
gradientsi^ Such modes may also couple to quantum 
systems^i"— . The emission in the vicinity of nanoanten- 
nas is also strongly modified with respect to free spaced!. 
Thus, to understand the dynamics of a quantum system 
it is essential to also account for these modified emission 
characteristics. As an example, a three-level system that 
is excited through an electric quadrupole transition and 
relaxes via two consecutive electric dipole transitions is 
investigated in detail in this contribution. 

The paper is organized as follows. At first, a mathe- 
matical framework to discuss the enhancement of higher- 
order multipole fields in the vicinity of optical nanoanten- 
nas is introduced. For the sake of definiteness, the discus- 
sion is exemplified at a suitably chosen nanoantenna, i.e. 
a nanoantenna which strongly enhances the quadrupole 
field. Next, the modification of excitation rates in quan- 
tum mechanical systems due to this quadrupolar en- 
hancement will be investigated. Finally, the hybrid sys- 
tem consisting of an optical nanoantenna and a quantum- 
mechanical three-level system will be studied in detail. A 
detailed appendix provides further explanations and ex- 
plicit derivations of formulas used in the main body of 
the manuscript. 

It is the purpose of this contribution to show that 
a properly designed nanoantenna can excite dipole for- 
bidden transitions in three-level systems due to the en- 
hanced higher-order multipole fields. It will be shown 
that the dynamics of the system are strongly altered by 
the presence of the nanoantenna and cannot be under- 
stood by the quadrupolar enhancement alone. Although, 
only shown here at the example of a three-level system, 
more complicated quantum systems can equally be con- 
sidered and the framework presented herein can be ad- 
justed to different possible experimental configurations. 



II. ENHANCEMENT OF HIGHER-ORDER 
MULTIPOLE FIELDS 

Previously, higher-order multipole transitions in 
hydrogen-like atoms have usually been assumed to be for- 
bidden (except for extreme situational^), since the cor- 
responding contributions as provided by an excitation 
field usually used in most cases, i.e. a plane wave, are 
much weaker than the dipolar ones. This can be esti- 
mated from a back of the envelope calculation that as- 
sumes some characteristic values. A stricter derivation 
can be found in, e.g., Ref. H^. The electric field of a 
plane wave varies spatially as exp [i(k ■ x)] w 1 -|- i(k • x) 
in the limit of k • x — 0. Here k is the wavevector that 
can be assumed as |k| w 10'' m"-'^ for visible light and the 
characteristic spatial extent of the atomic system being 
(|x|) « Zao w 10^^° m for hydrogen-like atoms. Only 
the first term in the expansion has been retained, cor- 
responding to the electric dipole component. In most 
cases this approximation is reasonable since the first or- 
der term in the Taylor expansion is three orders of magni- 
tude larger than the second term, which is attributed to 
both the electric quadrupole and magnetic dipole fields. 
Hence, for the given spectral domain and the usually con- 
sidered spatial extent of the atomic system, the excita- 
tion rates induced by the local quadrupole field of a plane 
wave are orders of magnitude smaller than those transi- 
tions induced by the electric dipole field. For this reason, 
quadrupole transitions are usually said to be inaccessible, 
i.e. they are forbidden. Clearly, components of octupolar 
or higher order are even weaker. 

In the presence of an optical nanoantenna the situation 
changes dramatically. Such plasmonic structures support 
highly localized near-fields that are characterized by huge 
gradients. Using a multipole expansion in spherical co- 
ordinates, the electric field in a coordinate system with 
origin Tq can be expressed as 

E (x, u) = y^}pmn (w; ro) N„„(x - tq, + 

m,n 

(w;ro)M„„(x - ro,w)] (1) 

with vector spherical harmonics Nm„ and Mm„ following 
the notation of Ref. [s^- Equation [T] is identical to a mul- 
tipole expansion in spherical coordinates (except some 
prefactors) with Pmn and qmn being the complex electric 
and magnetic multipole coefficients, respectively'^^. The 
order n — 1 corresponds to electric and magnetic dipoles, 
whereas n = 2 corresponds to quadrupoles and so on. 

For illustration purposes, in what follows the focus 
shall lie on the electric quadrupole. Moreover, to sim- 
plify the subsequent discussion, the more familiar Carte- 
sian quadrupole components Qij will be used. They can 
be obtained from the electric quadrupole coefficients pm2 
using linear transformations^^. 

As a referential optical nanoantenna that supports 
strong local electric quadrupole fields, an optical 
nanoantenna consisting of two strongly coupled silver 
nanospheres, sometimes termed as dimer, is investigated. 
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FIG. 2. (Color online) a) Integrated enhancement factor rjxz 
as a function of excitation wavelength for two separations 
of the dirtier nanospheres under plane wave illumination as 
sketched in b). The area Q pertinent for integration was 
chosen to be Az = 30 nm and Ax = 13 nm and fixed at 
y — [white box in b)]. A broadband enhancement of about 
six orders of magnitude for a separation of d = 3 nm can 
be observed which peaks at 370 nm and 437 nm. b) The 
local quadrupole enhancement ?7iT(ro) in the a::-z-plane for 
d = 3 nm and at A = 437 nm. c) Zoomed view of b) (color 
bar is maintained) where two symmetric stripes of maximum 
enhancement can be observed. 



enhancement factor rjij = r/l"'^ (tq) dV with respect to 
a certain domain n. Note that for enhancements obeying 
certain symmetries it is convenient to regard a lower di- 
mensional integration. Throughout this manuscript, the 
superscripts 'fs' and 'na' designate the free-space and the 
nanoantenna scenario, respectively. 

Figure [2] a) shows rj^z as a function of excitation wave- 
length for two different separations of the nanospheres. 
The dimer is illuminated by a plane wave as sketched in 
Fig. [21 b) . It can clearly be seen that for a separation of 
3 nm r]xz has a maximum of 1.6 x 10^ at 437 nm. Because 
it is strongest, we will exclusively focus on this resonance. 
For comparison Fig. [5] a) also shows rj^z for a nanosphere 
separation of 10 nm. The integrated enhancement fac- 
tor is approximately two orders of magnitude less than 
compared to the 3 nm separation. Therefore, it is ob- 
vious that the enhancement of the electric quadrupole 
contribution with respect to the near field stems from the 
strong coupling of the two nanospheres, which critically 
depends on their separation. All other components of r]ij 
are orders of magnitude less than rj^z and shall therefore 
be neglected here. Figure [2] b) shows 77^°"^ (ro) in the xz 
plane for a 3 nm separation at a wavelength of 437 nm. 
The enhancement of the electric quadrupole contribution 
takes place in a narrow spatial domain between the sur- 
faces of the two nanospheres. In Fig.[2]c), which displays 
a zoomed view of b) , it can be observed, that 1]^°^ consists 
of two symmetric stripes of approximately 5 nm width. 



The two silver nanospheres have a radius of 30 nm and 
are separated by either 3 nm or 10 nm. With such a sep- 
aration, quantum effects and/or a possible nonlocal ma- 
terial properties of the silver nanosphere may only con- 
stitute a minor contribution which we do not consider 
here^i^i; material parameters were taken from Ref. [H. 
The assumed geometry is in reach of state-of-the-art fab- 
rication techniques and can even be scaled to large ar- 
rays of strongly coupled nanospheresi^ It is well estab- 
lished that the scattering cross section of such a dimer 
structure exhibits a strong quadrupole contribution for a 
plane wave illumination direction parallel to the connect- 
ing line of the nanospheres'^^. However, the enhancement 
of local quadrupole fields, as probed by the quantum me- 
chanical system, was found to be strongest if the two 
nanospheres are illuminated perpendicular to the con- 
necting line (chosen to be the x-axis) with a polarization 
of the electric field parallel to it (see Figure. [5]). However, 
the quadrupole contribution to the far field is negligible 
for this illumination scenario. 

Considering plane wave illumination in the given coor- 
dinate system, the only nonvanishing component of the 
related quadrupole tensor in free space is Q^^^- This 
particular coefficient is linked to the multipole expan- 
sion coefficient as used in Eq. [T] via Q^z = — 
P12). Therefore, the local quadratic enhancement of a 
quadrupole field may be defined as the ratio 77]°'^ (ro) — 
I f 1 2 

\Qfj'{ro)/Q^z{ro)\ . One may also define the integrated 



III. FORBIDDEN TRANSITIONS 

After having considered the enhancement of the lo- 
cal electric quadrupole field near the nanoantenna, the 
resulting enhancement of quadrupole transition rates 
will now be examined. The transition rates from the 
i"^ to the j*"^ eigenstate of a non-relativistic quantum- 
mechanical system with unperturbed Hamiltonian Hq 
follows Fermi's golden rule 

r.. i^^i - ^) + '5 K + ^)} • (2) 

Here, the time-harmonic interaction potential Vi — 
g-icjtyj-^-j _^ e'"*y^(x) is assumed to be a small pertur- 
bation. Within this framework, the action of an electro- 
magnetic field on an electron with momentum operator 
p is given by the minimal coupling interaction potential 
Vi = —^A ■ p + h.c.^^, see also Appendix IB II Here, 
the weak-field approximation is applied and hence the 
ponderomotive potential shall be neglected. The 

electromagnetic fields are related to the vector potential 
A by B = curl A and E = iwA — gradC/. Choosing the 
Coulomb gauge where U = can always be achieved, 
the electric field is given by E = iwA. 

Furthermore, a decomposition of the interaction po- 
tential yields Vj = -i§ {[Hq, A ■ x] + [A, Hq] x}. In this 
form, the terms can be directly interpreted as electric 
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and magnetic coupling terms, V = Ve + Vm with 

K = -£[i?o,E-x],K^ = ^L.B (3) 

where fiB — ^ is Bohr's magneton and L = x x p 
denotes the angular momentum operator. A derivation 
of this decomposition can be found in Appendix IB 21 

Without loss of generality, the analysis shall be re- 
stricted to electric transitions. Referring to Eq.[T]and the 
relation Mmn(x) -x = 0, the electric interaction potential 
Ve can be fully characterized by the complex multipole 
coefficients p,„„ and the related vector spherical harmon- 
ics Nmn- Using the definition of the Nm„ following the 
notation of Ref. [30l one finds the explicit expression 



E • X = ^p,i 



n{n + 1) 



jn{kr)Pl^ (cos 6) exp(imi^) (4) 



where j„ denotes spherical Bessel functions and P™ as- 
sociated Legendre polynomials. The latter expression 
holds in the coordinate system of the quantum system, 
i.e. for To — 0. 

As mentioned earlier, this spherical representation can 
also be transformed to a Cartesian one and the mo- 
ments can locally be related to the coefficients of a Tay- 
lor expansion, E ■ x oc YJi...k Qi---kX^ ■ ■ ■ ■ Then, 
it becomes evident why the quadrupole coupling term 
(n = 2) is related to the near field gradients of the elec- 
tric field: Qij oc i9ii9j (E • x). However, care has to be 
taken regarding the normalization of the Qij which is not 
consistent throughout the literature. For the calculations 
in the main body of the manuscript, we have chosen the 
notation as given in Ref. |32| . 

Combining Eqs. [U - [31 the transition rates, related 
to different multipole orders, can be computed directly 
from a given field distribution. Hence, the transition rate 
at the atomic site ro, resulting from the electric poten- 
tial Ve, is proportional to the square modulus of a linear 
combination of the multipole coefficients Pmm i-e. 



r.,(ro) 



27re^ 



^Pmn(l"o) (i| N™„ • x|j) 



5 {uij ± u) 



(5) 

Equation [5] establishes the link between the enhance- 
ment of higher order electric transitions in a quantum- 
mechanical system and contributions of higher order mul- 
tipoles to the local field as investigated in Fig. [51 There- 
fore, the previously investigated optical nanoantenna can 
be used to increase the strength of electrical quadrupole 
transitions in a quantum mechanical system. 



IV. MODIFIED EMISSION 
CHARACTERISTICS 

After having demonstrated that near optical nanoan- 
tennas the quadrupole field can be enhanced by many or- 
ders of magnitude, the focus will now be directed to the 



effect of these multipole contributions on the dynamics 
of the entire quantum-mechanical system. For simplicity, 
a three-level system will be considered in the following, 
which is excited through an electric quadrupole transition 
and decays via two consecutive electric dipole transitions, 
see Fig. [31a). It must be noted that it is reasonably as- 
sumed herein that only the quadrupolar fields contribute 
to the quadrupole transition, i.e. the transition is dipole- 
forbidden. If this assumption does not hold, the inter- 
action to other enhanced multipole components of the 
field has to be taken into account following Eq. [5l As an 
example, a comparison to local dipolar enhancements is 
given in Appendix [XI In the subsequent discussion it will 
be shown how the emission characteristics of the entire 
system are modified in the presence of the nanoantenna. 
As a system of reference, a dimer with a distance of 3 nm 
between the silver nanospheres has been chosen. 



A. Dynamics of the Three-Level System 

The dynamics of the three-level system are governed 
by the following rate equations. They describe the pop- 
ulation of the three energy levels, i.e. 

= 710 • ni ~ • no , 
■hi = • n2 - 710 • ni , and 
"-2 = ro2 • no - 721 • n2 ■ 

Here, the denote the spontaneous decay rates from 
the i"^ to the j**^ level, whereas ro2 denotes the excita- 
tion rate of the quadrupole transition. The quadrupole 
transition takes place at a wavelength of A02 — 437 nm, 
for which the previously investigated optical nanoan- 
tenna provides the maximal enhancement of an electric 
quadrupole field [see Fig. [51a)]. 

In the system under investigation, quadrupolar emis- 
sion processes at A02 are neglected^, because the spon- 
taneous decay from the second to the first level is as- 
sumed to be much faster. The emission wavelength of the 
first dipole transition is assumed to be A21 = 3.47 /zm; 
consequently, the second dipole transition takes place at 
Aio = 500 nm. For both dipole transitions, stimulated 
processes are neglected. 

In free space and under plane wave illumination with 
unit intensity, the ratios rQ2/72i = 10~^ and 7{o/72i = 
10~^ were chosen. These rates correspond to real physi- 
cal systems as outlined in Appendix [Cl 

To investigate the enhancement of the efficiency of 
the quadrupole transition due to the nanoantenna, the 
amount of light spontaneously emitted at Aio ~ 500 nm 
per unit time shall serve as a figure of merit. There- 
fore, the rate equations are solved in equilibrium. In this 
regime, luminescence is given by h{^'^ = 7ig'^ni, where 
7j;g'^ is the spontaneous radiative decay rate. The total 
decay rate is given by the relation 710 = 7^0'^ -I- 7"™''^'^ 
where 75^^"''^'^ describes the nonradiative decay rate. 
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FIG. 3. a) Scheme of the quantum mechanical system that 
is placed in the vicinity of the nanoantenna. b) Local en- 

^ j_ r j_i 1 ■ • • rad,na / ■ rad,fs • ,^ 

hancement of the lummescence, i.e. /n^ , m the 

a;-z-plane. Because of saturation, the enhancement is not as 
strong as expected from the quadrupole enhancement alone, 
c) Luminescence enhancement as a function of the intensity 
at the location a; = nm, z — 4 nm. The dashed lines corre- 
spond to predictions for low and high intensities as discussed 
in section HV Bl fi) and (ii). 



Note, in this study in free space 7^g = Tio*^'^^, and 
^nonrad,fs _ q ^j^^ other hand, in the presence of the 

, nonrad.na • i j j_ ■ j t„ j_i -l 

nanoantenna 7^9 is merely determined by the ab- 

sorption in the metal. From the rate equations in steady 
state and 1 = uq + rii -\- n2, one can solve for ni finding 



■ rad 



721 Ti 



02 



710 721 + 721 + 710 To 



(6) 



which holds both in free space and in the vicinity of the 
nanoantenna. Since the quadrupole enhancement of the 
electric field varies spatially, the same holds for . Fur- 
thermore, the presence of the nanoantenna alters the lo- 
cal density of states. Thus, the spatial dependence of the 
spontaneous decay rates 7jj also has to be considered. 

Within the weak atom-field coupling regime, 
the spontaneous decay rates are given as 



where G is the 



component of the Green's tensor corresponding to 
the direction of the dipole moment dij of the specific 
transition^^. The Green's tensor is calculated at the 
position ro of the quantum system. Using this relation, 
the decay rates in the vicinity of the nanoantenna 7^"^ 
are related to the decay rates in free space 7^^ via 



7"77l- = S[G"^(ro,ro, 



[Gf^(ro,ro,c..,)] 



Since Qxz is the dominant quadrupole component, the 
emission from the dipole transitions is assumed to be 
either x- or z-polarized; hence is the arithmetic mean 
of these two contributions. 

In the subsequent discussion, ff{^'^ is chosen as the fig- 
ure of merit, assessing how efficiently the emission at 
Aio can be raised due to the enhancement of the elec- 
tric quadrupole fields. Figure [3] b) shows 



• rad.na / • rad,fs 



It can be seen that n^'^''^^ 



j • ra 



can be enhanced by over 
four orders of magnitude relative to free space. As ex- 
pected, in regions where the quadrupole enhancement is 
the strongest, the luminescence is strongly enhanced as 
well. However, due to saturation effects and nonradiative 
losses, characteristic values of the luminescence rates are 
about two orders of magnitude lower than the quadrupole 
enhancement. This can be understood in terms of differ- 
ent limiting cases of Eq. [5] and will be discussed in the 
following. 



B. The Effect of Excitation Intensities on the 
Luminescence Enhancement 



Naturally, the chosen ratios of the rates affect the lu- 
minescence enhancement. However, the excitation rate 
depends not only on the geometry but also on the 
intensity. Using Eq. [Sl one may estimate the lumines- 

1 , . rad.na / • rad.fs i , i • x •!„ 

cence enhancement /rii by the mam contribu- 

tions given in the nominator for free space and in the 
vicinity of the nanoantenna with respect to different in- 
tensities: 



• rad,na rad,na pna „,fs _i_ „, T^fs , -,fs pfs 

% ^ 7lO ^ £02 , 7lO 721 + 721 t 02 + 7lO J- 02 

■ rad,fs radjts "pfs _,na _i_ pna _i_ na pna ' v ' 

Ui 7io -1-02 7io 721 + 721 i 02 + 7io 02^ 

Purcell effect and dynamics of quantum system 

quadrupole enhancement 



In this form it tends to be evident that not only the en- 
hancement of the radiative rate of the luminescing tran- 

■ , • rad.na / rad,fs i ii 

sition, 7]^Q / 7]^Q , has to be taken into account as well 
as the quadrupole enhancement FqI /Fq2 but also the dy- 



namics of the quantum system given in the last term in 
Eq. [7l In passing we note that the enhancement of ra- 
diative rates in an environment is usually termed Purcell 
effect. 
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As discussed before, the excitation rate in free space 
at an intensity Iq was assumed to be {I = lo) — 
10~^72i- Also 721 = 10^7^0 1^^*^ been chosen implying 
^02 — ^o) ~ 10^^7io- Furthermore, in the vicinity 
of the nanoantenna 710 gets also hugely enhanced and 
710 ^ 721 holds. Equation [7] can be understood for dif- 
ferent limiting cases of low, intermediate, or high inten- 
sities. This is done in the following and will explain the 
results found in Fig. [3] in detail. 

(i) Low Intensities. First of all, one might consider 
the case of very low intensity, «C "fij both in free 
space and close to the nanoantenna. Then, one has 



rad,] 



rad,na 

7io 

• rad.fs rad.fs 

^1 7io ' 



pna 
^ 02 
pfs 
^ 02 



7^0 721 
7io 721 



racl.na 

Tio 

7?o" 



pna 
^ 02 

^ 02 



antenna efficiency x quadrupole enhancement 



since in free space 7^5 



^rad,fs^ One can see that the an- 
tenna efficiency 7io'^'"'^/7io is a limiting factor which will 
also naturally hold for the other limiting cases. Further- 
more, the luminescence enhancement is basically given by 
the quadrupole enhancement and can be calculated with- 
out consideration of the internal dynamics. Thus one can 
anticipate that the enhancement is most pronounced in 
this intensity regime. If the intensity is stronger, one can 
expect that the internal dynamics act as a bottleneck 
which will be examined in the following. 

(ii) High Intensities. Now the intensity will be as- 
sumed to be so high that in free space and in the vicinity 
of the nanoantenna S> 7y holds. Then, 
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consistent with Fig. [3] c) for I ^ Iq. Interestingly the 
enhancement is now independent of the quadrupole en- 
hancement. This results from the fact that in this inten- 
sity regime also in free space the excitation is the fastest 
process. 

(iii) Intensities comparable to Iq. Finally, one may 
consider an intermediate case where / « /q. Here, Fgj <^ 
721 but one may still assume that FqI 3> 721- This is the 
case in the central region of the nanoantenna which can 
be seen in Fig. [2] c) where the quadrupolar enhancement 
is shown to be in the order of 10^. Thus, 



. rad,: 



■ rad,fs 



rad,na 

7io 

rad.fs 

7io ' 



7?o72i 



pna 
^ 02 

pfs ^na p 
02 7lO ^ 



na 
02 



rad.na 

Tio 

7i"o' 



721 



antenna efiiciency x 10 



(8) 



It is important to note that the results outlined in 
Fig.[3]b) correspond to the discussed intermediate inten- 
sity limit. There, maximum luminescence enhancements 
in the order of 5 ... 6 • lO'^ were found. Looking at Eq. |S1 
this implies an antenna efficiency of approximately 50% 
in-between the nanospheres for the luminescence transi- 
tion. This could be confirmed by simulations shown in 
Fig.m 
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FIG. 4. The antenna efficiency 7io'*''^''/7io • It acts as a nat- 
ural limiting factor for the luminescence enhancement. Note- 
worthy, with respect to the actual form of the quadrupole 
enhancement outlined in Fig.[2]c), a perfect agreement of the 
luminescence enhancement predicted by equation Eq. [S] and 
the results given in Fig. [3] b) can be seen. 



It can be stated that under the given assumptions, the 
dynamics of the quantum system had to be taken into 
account to understand the luminescence enhancement of 
the system driven by a quadrupolar excitation. Only in 
the weak excitation limit, the luminescence enhancement 

T-G -L • rad.na / • rad.fs lor rad.na / na • 

simplmes to rij^ ' /^j^ = 7io /Tio - from tins 
relation it can clearly be seen that the local lumines- 
cence enhancement is directly proportional to the local 
quadrupole enhancement r]}^^. For a sufficiently strong 
excitation field, however, the entire process saturates and 
the overall enhancement decreases. 



V. CONCLUSION 

In conclusion, it was shown that optical nanoantennas 
can effectively enhance higher order multipole transitions 
which are typically considered forbidden in free space. 
This can be achieved by enhancing higher order multipole 
fields near the antenna. A quadrupole transition as the 
dominant excitation channel in a three-level system was 
considered. It was demonstrated how the enhancement 
of this transition can significantly intensify subsequent 
emission processes with respect to altered emission char- 
acteristics. Since the effects under consideration depend 
on geometrical parameters, the properties of the optical 
nanoantenna can be tailored and hence allow for direct 
implementations in spectroscopic schemes. It must, fur- 
thermore, be emphasized that the example of the isolated 
dimer is just a specific case of what is possible with ad- 
vanced optical nanoantennas. It can be anticipated that 
this work will give impetus for further research on such 
hybrid systems as well as plasmonic engineering. 
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APPENDIX 



Appendix A: Dipole vs. Quadrupole Enhancement 
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In the considered model system it was assumed that 
the energies of each transition are separated such that 
it is unlikely that different transitions are competing for 
a certain excitation. So, only the excitation rate with 
respect to the quantum system in free space was com- 
pared. One may also ask the question what happens if 
two transitions, namely a dipole and a quadrupole one, 
are energetically very close. Can one then make the lat- 
ter transition rate comparable or even stronger than the 
dipole transition? This should be possible when a cer- 
tain eigenmode of a structure can be tailored to have 
a strong quadrupole and only a negligible dipole field. 
For the dimer antenna under consideration this is not 
the case and was not intended. Here, also the electric 
dipole field gets enhanced. The result of a corresponding 
calculation is outlined in Fig. [51 There one can see that 
the electric dipole field gets enhanced approximately four 
orders of magnitude. This is still a few orders of mag- 
nitude smaller than the quadrupolar one. However, the 
difference for this particular antenna might not suffice to 
enhance a quadrupole transition to be stronger than a 
dipolar one. 

In general quantum systems, transitions might on the 
other hand be driven by a superposition of several multi- 
polar fields following Eq. [21 i.e. for Rydberg atoms. An 
enhanced quadrupolar field component then might allow 
an additional energy supply channel to such a system. 
However, as in the discussed three-level system, an en- 
hancement of higher order field componenents must be 
understood in terms of the whole dynamics of the quan- 
tum system. The approach might be most beneficial if 
this additional energy supply can account for an existing 
bottleneck for processes one wishes to enhance. 



FIG. 5. Local dipole enhancement factor (ro) = 

jd"'' (ro) (ro)| for the main dipole contribution dx- The 
configuration is the same dimer nanoantenna as used be- 
fore with the same excitation wavelength A = 437 nm. The 
dipole enhancement is approximately two orders of magnitude 
smaller than the quadrupolar one. 



straight forward, it will be outlined in the following two 
subsections. 



1. The Interaction Potential 

The action of an electromagnetic field on a given sys- 
tem within the framework of nonrelativistic quantum me- 
chanics is realized due to a replacement in the unper- 
turbed Hamiltonian in the following way: 

Hq (p, x) ^ (p - eA (t, x) , x) + eU {t, x) 
= Ho(p,x) + V/(p,x,t) 

which is called minimal coupling22ii2S . Throughout the 
analysis, a hydrogen-like Hamiltonian of the form 



Ho (p,x) = — p- 
2m 



1 Ze' 
47ren r 



(Bl) 



will be assumed. 

The relation of the fields to the electrodynamic poten- 
tials is given by 



Appendix B: Minimal Coupling 

In the main manuscript, the derivation of the mini- 
mal coupling potential Vj = — ^A • p was not presented, 
as well as its decomposition into an electric and mag- 
netic part given by Eq. [31 Although, Vi can be found in 
the literature and the derivation of its decomposition is 



B = curlA and E = -9* A - gradC/ . (B2) 
Choosing the Coulomb gauge in which 

divA and [/ = (B3) 
holds a p — and j = 0, one finds for the given Hamil- 
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H (p, X, t) = — (p - eA (x, t)f--. 

Am 47reo r 



1 
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[p^ — ep • A - eA • p 



+ eC/(x,t) 

1 Ze2 



where fiB = P = — i^V and L = x x p was used. So, 
it is natural to define 



m h e 
^^g^ ^ as the magnetic couphng. 



(B6) 



1 Ze' 



2rn 47reo J' 
+ ■ 



-Ap 

m 



Appendix C: On the Choice of Energy Ranges and 
(B4) Decay Rates 



where it was used that in the given gauge U — can 
always be achieved and p ■ A = A • p. Furthermore, the 
term in A^ has been neglected. This part corresponds to 
the PONDEROMOTIVE force. 



2. Electric and Magnetic Coupling - a 
Decomposition 

Now that the interaction potential Vi — —^A-p is 
known, it is desirable to further split it into an electric 
and magnetic part. This approach eases later interpreta- 
tion such as the attribution of transitions to electric and 
magnetic multipoles. 

First of all, one finds 

im r , 
A-p= —A - [iJo,x] 

ITTl 

= — {A ■ iJox - A • xi/o + HqA ■ X - HqA ■ x} 

n 

1 Ttl 

= — {[ffo,A.x] + [A,ifo]x} . 
n 

Now, one can state for the first term 

(to| [Hq, a ■ x] \n) = hujmn (™| A • X |n) . 

Because of E = — 0* A, this term can be interpreted in 
a time-harmonic dependency of the fields as the electric 
contribution. Then, E = iwA and the electric coupling 
may be introduced as 



Ve = 

On the other hand. 



e im , 
-[ffo,A-x] 

m n 

-^[i^cE.x]. 



(B5) 



1 



2m 



{Ap2 - p2A} 



A A = curlB 

2m 2m 

where the curl is acting only on B. It is further 

[A, Ho] X = curlB • x = (x x V) • B 

2m 2m 

= i^L.B 

e 



The numbers chosen for the calculations were moti- 
vated by order-of-magnitude estimations for quantum 
systems, namely hydrogen-like atoms and dye molecules. 
In this section, the choice of transition energies and used 
rates will be related to existing quantum systems. The 
suggested systems may be used in experiments. 



Energy Ranges 

A quadrupolar transition can be seen as the sum over 
all possible consecutive dipolar transitions. E.g. for the 
quadrupolar potential V (x) oc xy one finds 



Vij (x\{i\xy 



Hydrogen-like atomic systems shall be considered now. 
There, such transitions may be realized by the transition 
from an s to a d orbital using a p orbital as main inter- 
mediate step. These three levels then form an effective 
three level system as the studied one. 

The chosen energies correspond to such systems: 
Potassium has a quadrupolar transition X^a^M ~ 446 nm 
with intermediate step A4s^4p « 770 nm followed by 
Xip^M ~ 1.18 pm. The next alkali atom is Rubid- 
ium with X^s^id ~ 516 nm, X^s-^sp ~ 780 nm and 
X^p^id ~ 1.4 pm. For Cesium, one finds Aes->5rf ~ 
685 nm, Ags-s-ep ~ 894 nm and Agp-^sd ~ 2.9 pm. 

However, Lithium and Sodium have higher quadrupo- 
lar transition energies with wavelengths around 300 nm. 
Thus for such atoms an experimental realization may not 
be feasible - gold and silver lose their metallic character 
for such high photon energies. 



Decay Rates 

At a certain saturation intensity Is, dipolar excita- 
tion rates become comparable to spontaneous emission 
rates22ii£. For a two-level-system, Is is given by the in- 
tensity at which both states are equally likely. 

Hydrogen-like atoms exhibit an Is in the order of tens 
of W/m^ at optical frequencies. For instance. Rubidium 
has an Is = 16.4 W/m characterizing the 5s-5p transi- 
tion. On the other hand. Is may be more than 10^ W/m 
for dye molecules^. 
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At intensities above Is , the excitation rates for dipolar 
excitations become faster than their spontaneous coun- 
terparts. Quadrupolar excitation rates are for optical 
frequencies in the order of six to seven orders of mag- 
nitude smaller than dipolar ones. Thus, at intensities 
Iq w 10'^ "^/s, one might expect a quadrupole excitation 
rate to be in the order of (-^o) ~ 10~'^7io '^^^ 
sen in the calculations. 

The normalized intensity range in Fig. [3]c), I / Iq ~ 
10"'* . . . lO** corresponds to spectroscopic measurements. 
Intensities employed in pulsed systems are in the order 
of 10"'^^ . . . 10"'^'* W/m^ and for continuously operating sys- 
tems with plasmonic structures 10^ . . . 10^ W/m^. Hence, 



for dye molecules pulsed measurements are at an inten- 
sity of /pulse ~ 10® /s ~ 10^ ' '^/o and continuous wave 
measurements are at lew ~ /s ~ 10^"^ ^/q. As it was 
shown in Fig. [3] c), in-between these intensity ranges 
different limiting cases take place as discussed in sec- 
tion HVB] 

Furthermore, the fast nonradiative relaxation rate 
from state |2) to |1), was assumed to be the fastest 
process in free space. It was assumed to be independent 
of the environment, 'y2i = = lO^Tio- This assump- 
tion corresponds to fast thermal relaxations of excited 
atomic systems. 
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